Comparative study of quantum dynamics of a few bosons in a one-dimensional split 
hard-wall trap: exact results versus Bose-Hubbard-model approximations 
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We study the dynamical properties of a few bosons confined in an one-dimensional split hard 
wall trap with the interaction strength varying from the weakly to strongly repulsive regime. The 
system is initially prepared in one side of the double well by setting the barrier strength of the 
split trap to be infinity and then the barrier strength is suddenly changed to a finite value. Both 
exact diagonalization method and Bose-Hubbard model (BHM) approximation are used to study the 
dynamical evolution of the initial system. The exact results based on exact diagonaliztion verify the 
enhancement of correlated tunneling in the strongly interacting regime. Comparing results obtained 
by two different methods, we conclude that one-band BHM approximation can well describe the 
dynamics in the weakly interacting regime, but is not efficient to give quantitatively consistent 
results in the strongly interacting regime. Despite of the quantitative discrepancy, we validate 
that the form of correlated tunneling gives an important contribution to tunneling in the large 
interaction regime. To get a quantitative description for the dynamics of bosons in the strongly 
interacting regime, we find that a multi-band BHM approximation is necessary. 

PACS numbers: 34.50.-s, 31.15.ac, 03.75.Lm 



I. INTRODUCTION 

Bose-Einstein condensates in double-well potentials 
have attracted much attention in the past decades. As 
a paradigm model for studying the competition effect of 
quantum tunneling and interaction, the double well sys- 
tems have been widely studied from many aspects [l|- 
Ho| . Due to the experimental progress in manipulating 
ultracold atomic gases, both the trap potential and in- 
teraction between atoms can be implemented with un- 
precedented tunability ll|, and thus the dynamics of 



many-body quantum states of interacting bosons can be 
experimentally explored by loading the ultracold atoms 
in double wells. For the atomic double-well system, 
the atom-atom interactions in Bose-Einstein condensates 
have been found to play important roles in the dynam- 
ics of the system. The competition between tunneling 
and interaction leads to many rich and interesting ef- 
fects, such as the Josephson oscillation and self-trapping 
phenomena [H-0|- Moreover, novel correlated tunneling 
dynamics in interacting atomic systems characterized by 
a small number of particles has recently been observed 
experimentally (l2j , which attracted particular attention 
in the study of the dynamics of few-atom systems fl3l . [l4j ]. 

So far most of theoretical works for double-well sys- 
tems are based on the two- mode approximation l6|j . 
For atoms confined in a double well potential, if all of 
them are prepared in one well initially, they will oscil- 
late in the form of Josephson oscillation for weak enough 
interaction, or they will stay in one well (self-trapping) 
when the interaction is above a critical value. Both the 
Josephson oscillation and self-trapping phenomena can 



be understood within the two-mode approximation and 
have been observed experimentally in cold atomic sys- 
tems [16]. However, if interactions are strong enough, 
the mean-field theory and the two-mode approximation 
are not expected to be valid as higher orbits are occu- 
pied. In the strongly interacting regime, strong interac- 
tions between atoms may fundamentally alter the tunnel 
dynamics and result in a correlated tunneling, which was 
explored most recently in ultracold atoms [l3|, [l4| • The- 
oretically, the correlated paring tunneling was studied 
by multi-configuration time-dependent Hartree method 
[l3l [l5| and also in the scheme of the extended Bose- 
Hubbard model with an additional term of correlated pair 
tunneling [141 ]. In order to understand the dynamics of 
double well systems from weakly to strongly interacting 
regime in an unified scheme, in this work we study the 
dynamical properties of a few bosons confined in a one- 
dimensional (ID) split hard wall trap with the repulsive 
interaction strength varying from zero to infinity. Exper- 
imentally, the effective interaction strength can be tuned 
by using Feshbach resonance or the confinement-induced 
resonance to the strongly interacting Tonks-Girardeau 
(TG) limit pj], [3, which makes it possible to explore 
the novel dynamics even in the TG limit. 

In the TG limit, the bosonic systems exhibit the fea- 
ture of fermionization [H3H24T ] . In the strongly inter- 
acting regime, the mean field theory generally fails to 
describe the properties of fermionization. In order to 
characterize the crossover from weakly interacting con- 
densation to strongly interacting TG gas, some sophis- 
tical theoretical methods, such as the exact diagonaliza- 
tion method [H, [21 
Hartree method 
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and multi-orbital self-consistent 
231 124| have been applied to study the 
static few-boson systems. In this work, we shall apply 
the exact diagonalization method to study the dynam- 
ical problem in the ID double-well system. The exact 
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diagonalization method can produce numerically exact 
results and allows us to give an unified description for 
both the weakly and strongly interacting regime. As a 
comparison, we also investigate the dynamics based on 
the two-site Bose-Hubbard model by considering both 
the two-mode and muti-mode approximations. Compar- 
ing the results obtained from different methods, we con- 
clude that one-band (two-mode) BHM approximation is 
efficient to describe dynamics in small interaction regime, 
but the multi-band BHM approximation is needed if we 
want to describe the dynamics of bosons with large in- 
teraction quantitatively. We also validate that the form 
of pair tunneling gives an important contribution to tun- 
neling in the large interaction regime. 



II. MODEL AND METHOD 

We consider a few bosons with mass m confined in an 
one-dimensional split hard wall trap, which is described 
by the Hamiltonian (K = m = 1) 



H = / ${x) 



1 d 2 

w + v{x) 



k8(x) 



ip(x)da 



c / ■)/>' (x)ip* (x')5(x — x')i/j(x)i/)(x')dxdx' . (1) 



Here V(x) is a hard wall trap which is zero in the region 
(—a, a) and infinite outside, k is a tunable parameter 
which describes the strength of zero-ranged barrier at 
the center of the trap, and c is the interaction strength 
between particles determined by the effective ID s-wave 
scattering length. Here the double well is modeled by 
the ID split hard-wall trap with a 5-type barrier located 
at the origin and the tunneling amplitude between the 
left and right wells can be tuned by the barrier strength 
k [H, [l(|. To study the tunneling dynamics, the barrier 
strength k is initially set to be infinity and the system 
is prepared in the ground state of the left well. At time 
t = 0, we suddenly change k to a finite value and study 
the dynamical evolution of the initially prepared system. 

For c = 0, the single particle stationary Schrodinger 
equation associated with the Hamiltonian ([1]) can be 
written as 



1^ + V(x) + k6(x) 



ip n (x) = e n ip n (x), (2) 



where (p n (x) are the complete set of orthonormal eigen- 
f unctions and e„ the corresponding eigenenergies. Here 
n = 1, 2, 3, • • • gives the ordering number of the single- 
particle energies. According to the parity symmetry 
of the eigenfunctions, the state tp n (x) is symmetric for 
odd n (n = 1% — 1) and antisymmetric for even n 
(n = 2i). The single-particle energies are ordered alter- 
natively corresponding to symmetric and antisymmetric 



states. For convenience, we also represent ip2i-i(x) = 
tfi } s{x) and ifi2i(x) — ifii t A(x) with the subscript S 
(A) indicating the symmetric (antisymmetric) func- 
tion. The single-particle antisymmetric eigenfunctions 
are <p^a{x) = -^sin(^), i — 1,2,3,... with their cor- 
responding eigenenergies e^A = = [iit /a) 2 /2, and the 
single-particle symmetric eigenfunctions are tfi_s(x) — 
C[cos(px) - ~ sm(px)]8(~x) + C[cos(px) + % sm(px)]9(x) 
with their corresponding eigenenergies e^s = £2i-i = 
p 2 /2, where the wave vector p is determined by tran- 
scend equation p/ n + tan(pa) = and C is the normal- 
ization constant. It is true that the barrier only influ- 
ences symmetric eigenfunctions, but does not influence 
antisymmetric eigenfunctions for any barrier strength k. 
Further, the (2i — l)-th eigenenergy is close to the 2i- 
th eigenenergy (i — 1, 2, 3...) gradually with the increase 
of barrier strength and they become degenerate in the 
limit k — y oo . For simplicity, we set a = 1 and discuss 
a large barrier strength k — 50 in this paper. In this 
case, single-particle eigenenergies in split hard wall trap 
are ei = 4.74341, e 2 = ir 2 /2, e 3 = 18.9764, e 4 = 2tt 2 , 
£5 = 42,7073,... respectively, see FigflTa). In contrast 
to the small energy gap between the (2i — l)-th state 
and the 2z-th state, there is a relatively very large energy 
gap between the 2i-th state and (2i + l)-th state, i.e., 
£i,A — &i,S "C £i+l,S " £i,A- 

Expanding field operators as *S>(x) = XmLi L Pn(x)a n , 
the many-body Hamiltonian ([1]) takes the form 



H = t n a) n a n +c ^ InmpqO^a^apaq, (3) 



where at(a„) is bosonic creation (annihilation) opera- 
tor for a particle in the single particle energy eigenstate 
(f n . The interaction integral parameters l nmvq are cal- 
culated through I nmp q = Jl a Vn{x)¥m{x)ip p {x)ipq{x)dx. 

The eigenstate of this Hamiltonian can be obtained by 
numerical exact diagonalization in the subspace of the en- 
ergetically lowest eigen-states of a noninteracting many- 
particle system (Ioli2l|. 

When the barrier strength k is large, the split hard 
wall trap can be considered as a double well. Similar to 
the case of optical lattices, the local Wannier functions 
Wl(x) {Wft(x)) at the left (right) well with the energy 
band indices i can be defined as W l L {x) — l/v / 2((p ij 5(x) + 
ip it A(x)) and W l R {x) = l/y/2((p itS (x) - ^ A {x)). From 
the symmetry of tpi^s{x) and ipi t A(x), one observes that 
Wl(x) = W l R {— x). If we expand the bosonic field oper- 
ator as, tp(x) = J2i a i,LWl(x) +J2i a i,R w h( x ) > where 
a i,L(R) is the bosonic annihilation operator for a particle 
at left (right) well, the Hamiltonian can be written as the 
form of two-site Bose Hubbard model 
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i,j i,j i,j,k,l u,fi,~t,5 



where the integral J", = da;(W^(a;))*i?oW|(x), 

with i?o = ~~2'§x I + and the interaction in- 

tegral [T'J^ = c J dx(W^x))* (W^x)yW^(x)W^x). 
The subscripts a, (9,7,5 £ are the well indices 

and the superscripts k, lG {1, 2, 3, ...} are the energy 
band indices. Here we note J RR = J^ L for the symmet- 



ric double well. The Haniiltonian Q can be divided into 
intraband and interband parts, that is 

H = Hi + Hinterband- (5) 

i 

The i-th intraband Hamiltonian can be written as 



Hi = (e hL nix + e ijR n itR ) + [Ji + 2{n hL + n hR - 1) J-] (at L a itR + a\ R a i}L ) + 

C^oK,i( n »,i _ 1) + n i,R,{ n iM - 1)] + U l LR (a\ L a\ L a itR a i:R + aj R a\ R a it La lt L + 4:Ti ltL n ltR ), (6) 

I 



where n i:L = a] L a hL , n iiR = a\ R a iiR , e i:L = J% L , e iiR = 

Jrr> Ji = J % lr = Jrl IS ^e intraband hopping energy 
between left and right wells, J[ — U"^ LR — U RRRL , 
= ^llll = U RRRR is the on site interaction energy 



and U % LR = U % l % l RR is the intraband pair hopping energy. 
It is easy to check that e^j, = 6i <R — (e^s + e, l ^)/2 = ^ 
and Ji = (e^s ~ £%,a) /2. The interband Hamiltonian 
reads 



Hinterband = ^ E t/ ^f/3,7,«( J,a a i, i 8 a i.'Y i,« + a z> a i/3 a J>7 a ^) + E XI U ^y,S a i, a a ],p a k,yai,S (?) 
a,/3,7,<5 i,j,k,l a,/3,7,5 



where C£? ft7j4 



and the summa- 
tion in the second part contains three-band terms 
where only three different energy band indices exist and 
four-band terms with i / j / ); ^ 1. Most inter- 
band interaction terms are very small except for terms 
of U^ LLL and U RRRR defined on the same site and be- 
tween the ith and jth bands. For the split system with 
barrier strength n — 50, we have c^l = £i,R = 4.83911, 
Ji = -0.0957, Ul = 1.48446c, J{ = -0.00366c, U X LR = 
0.0003c and U^lll = 0.98866c, while the interaction 
strengths containing three or four energy band indices 
are very small, for example Uj^l L = —0.000969c and 
'/'f /;•/,■ 0.00107c. 

Now we turn to consider the dynamical behavior for 
the iV-boson system. Initially (for t < 0), the barrier 
strength k is set to infinity and TV bosons are prepared 
in the ground state of the left well. In this case, the 
split hard-wall trap actually reduces to two separated 
hard-wall traps, so the eigenvalues and eigenfunctions 
of a single particle in the left well are (7rn/a) 2 /2 and 
fn(x) = —yj2/a sin(n7re/a) , with n — 1,2,3... respec- 



tively. Then at t — 0, K is suddenly changed to a finite 
value, for example k — 50 in the present work, and the 
corresponding Hamiltonian is Hf = H(k = 50). After k 
is changed, the time-dependant wave function is given by 



\9(t)) =e-^ t |*(0)) 



OO 

E 



C n e- iE ^ r , 



(8) 



in which weight coefficients are C n — ($ n | v I'(0)) and they 
satisfy normalization condition J^Li = 1- Here $„ 
and E n are the eigenstates and eigenvalues of Hf, respec- 
tively. In order to see how the initial state trapped in the 
left trap evolves, we shall use revival probability 



F(t) = |(*(i)|*(0))| : 



E 



C 2 n C 2 m sm 2 [(E n -E m )t/2}, (9) 



the reduced single- particle density matrix 

p(x,x',t) = {*(t)\&(x)*(x')\*(t)), (10) 
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FIG. 1: (a)Single-particle energy levels for the split hard wall 
trap with k = 50. (b) The eigenenergy of two interacting 
bosons changing with the interaction strength c in the split 
hard wall trap with k = 50. 



and pair correlation function 

g (2 \x u x 2 ,t) = (^(t)\^(x 1 )^(x 2 )^(x 1 )^(x 2 )\^(t)) 

after time t to describe the dynamics in split hard-wall 
trap system. 



III. RESULTS AND DISCUSSIONS 

Before studying the quantum dynamics of many-body 
systems, we first recall the tunneling dynamics of a single 
atom. If there is only one boson in this split hard wall 
trap, the initial state is just the ground state of left well, 
that is 'F(O) = — V2 silica;) with energy 7r 2 /2. After 
the barrier strength k switches on to a finite but large 
strength, the weight coefficients of ground and the first 
excited state of the finial Hamiltonian Hf are C\ ss \/2/2 
and C2 ~ \/2/2. At time t, the wavef unction reads 



= e-^ t [cos(J 1 t)\Wl) +isin(Jit)|WA)],(12) 

where /ii = (ei,s + ei,yt)/2. It is obvious that the bo- 
son stays in the left well with the probability of cos 2 (Jii) 
whereas in the right well with the probability of sin 2 (Jit). 
Consequently, the boson oscillates back and forth be- 
tween two wells with period r = 2ir/(ei t A ~ e i,s) = 
— 7r/Ji, which is influenced by the barrier strength k 
through controlling the energy gap between ground state 



and the first excited state. Correspondingly, the fidelity 
F(t) ps cos 2 (Jit) oscillates periodically between 1 and 0. 

For a many-body system, no an analytical expression 
like Eq. (fT2|) is available. Nevertheless, when the atom 
number is small, we can resort to the full exact di- 
agonalization method to calculate the energy spectrum 
and eigenstates via directly diagonalizing the Hamilto- 
nian ([3]). Consequently the time-dependent wavefunc- 
tion, revival probability, single-particle density matrix 
and correlation function are straightforward to be cal- 
culated via Eq.([8|- (1TT1) . For a continuum system, we 
need truncate the set of single-particle basis functions 
to the lowest L orbitals (modes) and the basis dimen- 
sion of a TV-particle system with L orbitals (modes) is 
given by D = (N + L - 1)\/[N\ (L - 1)!]. In general, 
one needs L ^> N and it is a formidable task to get the 
full spectrum of the many-particle system as the parti- 
cle number N becomes large. Therefore, despite the fact 
that the exact diagonalization method can be applied 
to deal with the interacting boson systems in a numeri- 
cally exact way for all relevant interaction strengths, it 
only restricts to small particle systems. When the in- 
teraction strength is weak, the two-site Bose-Hubbard 
Hamiltonian under single-band (two-mode) approxima- 
tion is widely taken to be the model system for the study 
of the dynamics of the double-well system. One of ad- 
vantages of the two-site Bose-Hubbard Hamiltonian Q is 
that every term in the Hamiltonian has a straightforward 
physical meaning which can help us to understand the 
physical consequence of different terms. Furthermore, 
in the scheme of the two-site Bose-Hubbard model, the 
system is much more tractable both analytically and nu- 
merically and a large system can be studied. In the fol- 
lowing, we shall first present exact numerical results by 
exact diagonalization and then results based on the two- 
site Bose-Hubbard model under two-mode (single-band) 
and multi-mode (multi-band) approximations. 



A. Exact result by exact diagonalization 

We first consider the two-boson case. If two bosons 
are initially prepared in the ground state of the left well 
as the initial state of the system, which can be gotten 
by exact diagonalization method. Through diagonaliz- 
ing the second quantized initial Hamiltonian Hi n in the 
Hilbert space spanned by the single-particle eigenstates, 
we get the initial state ^(0). Similarly the eigenenergy 
and eigenvectors for Hf can also be obtained. As an 
example, we plot the lowest five eigenenergies of two in- 
teracting bosons in the split hard wall with k = 50 versus 
the interaction strength c in Fig. [ljb). 

When the interaction is absent, bosons just oscillate 
back and forth between two wells and return to their ini- 
tial state after a Rabi period r, which is the same as 
one boson's. For interacting bosons, there will be many 
differences. The revival probability F(t) as a function 
of time t is shown in Fig|5] for various c. When inter- 
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action strength c is very small, the revival probability 
F(t) still displays the oscillating feature and the system 
returns to their initial state with probability close to 1 
after a longer period. The revival time becomes longer 
with the increase of the interaction strength c. When c 
reaches a certain value (c ~ 5), F(t) approaches 1 with 
tiny oscillations within a very large time scale, which is 
known as the self trapping phenomena. In this regime, 
the tunneling to the right well is dynamically suppressed 
and two bosons stay in the left well stably. As the inter- 
action strength c increases further to the stronger regime, 
F(t) begins to decrease more quickly and two bosons can 
tunnel to the right well again. In the limit of c — > oo, 
F(t) approaches zero quickly and then oscillates between 
and 0.521, finally it approaches 1 after almost a Rabi 
period. To see clearly how the atoms tunnel between 
the left and right wells, we display the corresponding 
time-dependent reduced single particle density matrix 
p(x,x',t) in FigJ3] for several typical c. The diagonal 
contribution p(x, x' , t) along x — x' is just the single par- 
ticle density distribution. In FigUl the pair correlation 
function g^ 2 \x\, X2, t) are also displayed. The pair cor- 
relation function g^(xi,X2,t) shows the probability of 
finding one particle at point x\ and another particle at 
point X2 in one measurement. As shown in (al), (a2) and 
(a3) of Figj3] and FigfJ] the non-interacting bosons can 
tunnel from left well to right well, and then go back to left 
well, which forms a period of Rabi oscillation. However, 
as shown in (61)-(63) of FigJ3]and FigfJ] the single parti- 



FIG. 3: (Color online) Reduced single-particle density matrix 
p(x,x' ,t) of two interacting bosons for (al-a3) c = 0, t = 
0,5,15, (bl-b3) c = 5, t = 0,5,15, and (cl-c6) c = oo,t = 
0, 5, 16, 21, 22.5, 32. Each plot spans the range — 1 < x, x' < 1. 



cle density distribution and the pair correlation function 
g( 2 \x, x', t) have no obvious change within a Rabi pe- 
riod, and no oscillation between the left and right traps 
is observed in this self trapping regime with c = 5. While 
in the fermionization limit, the oscillation phenomenon 
appears again in (cl) — (c6). Both of two bosons tunnel 
to right well when t = 21 (see (c4)) and go back to the 
left well for t = 32 (see (c6)). Our results are consis- 
tent with results in [131 based on the multi-configuration 
time-dependent Hartree method. 

Next we consider the dynamics of systems with more 
bosons. The dynamics of N = 3 system is similar to the 
N = 2 case, except that the system enters self-trapping 
regime earlier than N = 2. As shown in Fig. [5l when 
c = 0.5, the system already displays the feature of self 
trapping with the fidelity F(t) ~ 1 with tiny oscillations 
within a very large time scale. Similarly, in the TG limit, 
bosons can tunnel to the right well more easy and return 
to the left well approximately after a Rabi period. The 
tunneling dynamics of hard-core bosons is very similar to 
their correspondence of free fermions (28|. In FigjSl we 
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FIG. 4: (Color online) The pair correlation function 
g (2) (xi,a;2,t) for (al-a3) c = 0,t = 0,5, 15, (bl-b3) c = 5,t = 
0, 5, 15, and (cl-c6) c = oo, t = 0, 5, 16, 21, 22.5, 32. Each plot 
spans the range — 1 < x,x' < 1. 



FIG. 6: The revival probability F(t) of the JV-boson system 
changes with time t with infinity interaction. 

plot the revival probability F(t) for systems with N — 
1 to N = 4 in the TG limit. It is shown that there 
is an obvious peak around the Rabi oscillation periods 
for various N, which implies that the system can return 
to the left well with probability close to 1 after a Rabi 
period. One can understand this from the Bose-Fermi 
mapping, i.e., bosons in the TG limit can be mapped 
into a spinless free Fermi system [28]. For the case with 
k = 50, we can check that Ji w i 2 Ji for i — 1,2,3,4. 
The N atoms initially occupy the iV-lowest single-particle 
levels of the left well, and roughly speaking, each particle 
tunnels with the Rabi period Tj = —ir/Ji. However, when 
the particle number N becomes large, the relations of 
Ji fa i 2 J\ break down and even the dynamics in the TG 
limit can be quite complex. 




FIG. 5: The revival probability F(t) of three particles system 
changes with time t with different interaction. 



B. BHM approximation 

Now we consider the two-site BHM described by the 
Hamiltonian (U). If the interaction strength is much 
smaller than the level spacing between the first band 
and the second band defined as A = ^2 — fJL%, one may 
expect that the system can be approximately described 
by the single-band BHM. Under the one-band approxi- 
mation, the Bose- Hubbard model is described by Eq.© 
with i = 1. We note that the pair hopping term U\ R in 
Eq. © is generally very small in comparison with the on- 
site interaction, for example, in the present work we have 
J7q f=a 4948L r ^ i? . Therefore in many previous works, the 
pair hopping term is omitted and a simplified single-band 
BHM given by 

H = (ii(ni tL + ni >R ) + Ji(a\ L ai jR + a\ M a ltL ) + 

Uo - 1) + - 1)], (13) 

has been widely used [j. \2^2l\ . Here Ji = Ji + 2(n ltL + 
ni.R — Since U\ R <C J7g in the whole interacting 
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FIG. 8: The revival probability F(t) for the two-boson system 
changes with time t under one-band BHM approximation with 
(a) Ul R = and J[ = 0. (b) U\ R = 0. 



FIG. 7: The revival probability F(t) for the two-boson system 
changes with time t under one-band BHM approximation. 



regime, the pairing hopping term is not expected to sig- 
nificantly change the static properties. However, when 
the term U\ R is comparable with the hopping amplitude 
Ji, it may give significant contribution to the dynamics, 
which has been emphasized in Ref . [lij ■ In the weakly in- 
teracting regime, the term of J[ is also usually neglected 
as its revision to hopping energy can be attributed to J\ . 
As we shall illustrate later, when the interaction strength 
is very large, the contribution of J[ can not be neglected 
since J{ cx c. 

By using the single-band BHM approximation, we 
study the tunneling dynamics of the two-site BHM de- 
scribed by Hi given by Eq.© with all particles prepared 
in one site initially. The revival probability F(t) chang- 
ing with time t under one-band BHM approximation are 
shown in FigfT] for the two-boson system. As shown in 
FigLZK, when the interaction strength is not very large, 
the single-band BHM gives quantitatively consistent de- 
scription of the dynamics in comparison with the exact 
results by exact diagonalization (see Figf^k). With fur- 
ther increase in the interaction, as shown in Figj7]3, al- 
though the single-band BHM with pairing hopping term 
can describe correctly the enhancement of tunneling, it 
does not provide quantitatively consistent results in com- 
parison with results of exact diagonalization in Figj2j). In 
order to see clearly the effect of the pair-tunneling term, 
we also study the dynamics governed by the simplified 
BHM of Eq. (IT51) without the pair-tunneling term. To see 



the effect of the term of J[, we consider both cases for 
the Hamiltonian of Eq. (|13[) with or without the term of 
J[. When the interaction strength is not so strong (for 
example, c < 5), we find that the results are almost the 
same as that presented in Fi^T^. That means that the 
terms of pair tunneling and J[ are not important when 
the interaction is weak. However as shown in Fig|8] (a) 
and (b), the dynamics in the strongly interacting regime 
shows quite different behaviors if the pair tunneling term 
is absent. Comparing Fig|7jD and FiglHl we can conclude 
that the pair tunneling term of V\ R gives an important 
contribution to tunneling in the large interaction regime. 
Comparing FigfSJi and FigfSp, we find that the term of J[ 
also plays an important role in enhancing the tunneling. 

The dynamics for the three-boson system is shown in 
FiglHl Comparing with the exact dynamical results in 
Figf51 we find that the one-band BHM approximation 
can describe the dynamics well only when the interac- 
tion strength is small so that Uq <C A. In contrast to 
the two-boson system, the pair-tunneling term has less 
significant effect on the enhancement of the tunneling. 
For very large c, although the system can tunnel to the 
right well, it does not give quantitatively consistent re- 
sults in comparison with results of the numerical exact 
diagonalization. 

From the above results, we know that one-band BHM 
approximation is not enough to give a quantitatively de- 
scription for the dynamics of interacting bosons in the 
large interaction regime. To get better results, we need 
keep more band levels and use the multi-band BHM given 
by Eq.(|l]) in our calculation. In FigfTUl we show the 
ground energy of the two bosons versus the interaction 
strength using exact diagonalization method, one-band, 
two-band, three-band and four-band BHM approxima- 
tions, respectively. It is shown that one-band BHM ap- 
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FIG. 9: The revival probability F(t) for the three-boson sys- 
tem changes with time t under one-band BHM approximation. 

proximation can describe the ground energy very well 
when the interaction c < 1, and two-band BHM approx- 
imation can describe well in the region c < 10, while 
three-band and four-band BHM approximations are ef- 
ficient to describe the ground energy of two bosons well 
even for c = 100. For the dynamics problem, in order to 
get a quantitatively consistent results with the exact di- 
agonaliztion results, we find that more bands are needed 
in comparison with the static problem. In Fig llll we 
display the results of F(t) for the two-particle systems 
with various c within the multi-band BHM approxima- 
tion. As shown in the figure, the result based on a five- 
band BHM approximation for c = 10 already quantita- 
tively agrees with the exact numerical result. For c = 50, 
a ten-band BHM approximation is required for a quan- 
titatively consistent result. The result for c = 300 based 
on an eighteen-band BHM approximation is also given 
in FigOT] In comparison with Fig. 2b and Fig. 7b, we 
find that there exists only a qualitative agreement with 
the exact diagonaliztion result although it is much better 
than the result of the single-band BHM approximation. 

IV. SUMMARY 

In summary, we have studied the dynamical proper- 
ties of a few bosons confined in an one-dimensional split 
hard wall trap by both the exact diagonaliztion method 
and the approximate method based on the two-site Bosc- 
Hubbard model. The system is initially prepared in the 
left well of the trap by setting the barrier strength of 
the split hard wall trap to infinity, and then it is sud- 
denly changed to a finite value. With the increase in the 
interaction strength of bosons, the system displays the 



Josephson-like oscillations, self trapping and correlated 
tunneling in turn. Comparing results obtained by two 
different methods, we conclude that the one-band BHM 
approximation can quantitatively describe the dynam- 
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FIG. 10: The ground energy for the two-boson system ob- 
tained by exact diagonalization method and i-band BHM ap- 
proximation with i = l,2, 3, 4. 




FIG. 11: The revival probability F (t) changes with time t un- 
der multi-band BHM approximations, including a five-band 
BHM approximation for c = 10, a seven-band BHM approxi- 
mation for c = 20, a ten-band BHM approximation for c = 50 
and a eignteen-band BHM approximation for c = 300. 



ics in the weakly interacting regime, but the multi-band 
BHM approximation is needed if we want to describe the 
dynamics of bosons with large interaction quantitatively. 
We also validate that the form of correlated tunneling 
gives an important contribution to the tunneling dynam- 
ics in the large interaction regime. 
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